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Using first-principles tecliniques, we calculate the renormalization of the electron Fermi velocity 
and the vibrational lifetimes arising from electron-phonon interactions in doped bilayer graphene 
and in graphite and compare the results with the corresponding quantities in graphene. For similar 
levels of doping, the Fermi velocity renormalization in bilayer graphene and in graphite is found to 
be approximately 30% larger than that in graphene. In the case of bilayer graphene, this difference 
is shown to arise from the interlayer interaction. We discuss our findings in the light of recent 
photoemission and Raman spectroscopy experiments. 



Since the fabrication of crystalline graphitic films with 
a thickness of only a few atomsi^i^i^i^, single- and 
double-layer graphene have received considerable atten- 
tion-. These materials are promising candidates for na- 
noelectronics applications because of the high mobility 
of charge carriers in these systems and the tunability of 
their electronic properties by gating^. Since electron- 
phonon (e-ph) interaction plays an important role in the 
dynamics of charge carriers^'*, understanding its effects 
in single- and double-layer graphene is of crucial impor- 
tance for graphene-based electronics. 

The e-ph interaction in metals modifies the dynamics 
of electrons with energy near the Fermi level by increasing 
their mass and reducing their lifetime. The mass renor- 
malization can be described in terms of the e-ph coupling 
strength A„k, defined as the energy derivative of the real 
part of the phonon-induced electronic self-energy TinkiE) 
at the Fermi level Ep: A„k = — 9ReI]„k(-E')/9i?|_E=_EFJ 
where n, k and E are the band index, the wavevector 
and the energy of the electron, respectively^. The elec- 
tron mass renormalization can be obtained from the e-ph 
coupling strength through m* /m =1-1- A„k where m and 
m* are the bare band mass and the renormalized mass, 
respectively. The e-ph interaction also gives rise to a 
phonon lifetime Ti,q = h/T^^, where Fj^q is the phonon 
linewidth, i. e. , twice the imaginary part of the phonon 
self-energy arising from the e-ph interaction. Here v, q 
and u are the phonon branch index, the wavevector and 
the energy of the phonon, respectively^. 

These quantities can be calculated from first-principles 
within the Migdal approximation as^ 



A„k = V' / l5m«,t/(k,q) 
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Here e„k and ujq^ are the energy eigenvalue of an elec- 



tron with band index n and wavevector k and that of a 
phonon with branch index v and wavevector q, respec- 
tively. Abz is the area of the first Brillouin zone where 
the integration is performed. The quantities /„k and 
TZjyq are the Fermi-Dirac and the Bose-Einstein factors, 
respectively, and g„„,,,(k, q) = (mk + q| AKq(r) |rtk) 
is the scattering amplitude of an electronic state |nk) 
into another state |mk + q) resulting from the change in 
the self-consistent field potential AViyqi^) arising from a 
phonon with the branch index v and the wavevector q. 

Electron wavefunctions and energy eigenvalues are ob- 
tained using ab initio pseudopotential density functional 
theory calculations^'' within the local density approxi- 
matioiiiiii^. Phonon frequencies and eigenstates are ob- 
tained through density functional perturbation theory^. 
We have used a planewave basis with a kinetic en- 
ergy cutoff of 60 Ry. The core- valence interaction is han- 
dled using norm-conserving pseudopotential a^^'^^ . The 
integration in Eq. ([1]) for graphene and bilayer graphene 
is performed by summation over 300x300 points and for 
graphite it is performed with 90x90x30 points in the 
irreducible part of the Brillouin zone. The integration 
of Eq. (m for graphene and bilayer graphene is done by 
summation over 1000 x 1000 points in the irreducible part 
of the Brillouin zone, and the 6 function is replaced by 
a Lorentzian with 15 meV broadening for convergence. 
Electron and phonon wavefunctions, energy eigenvalues 
and the e-ph coupling matrix elements in these extremely 
dense grid sets are obtained by a recently developed inter- 
polation schemeiii^ based on maximally localized Wan- 
nier functionsi^i^. The Fermi-Dirac and Bose-Einstein 
factors are evaluated at the temperature T = 15 K in all 
the calculations. Charge doping is modelled by adding 
or removing electrons from the simulation cell and by 
using a neutralizing background. In this work we as- 
sume that the layers in bilayer graphene and in graphite 
are arranged according to the Bernal stacking sequence 
[Fig.Ha)^. 

Figure [T] shows the e-ph coupling strengths A„k in 
graphene, bilayer graphene and in graphite calculated 
along the reciprocal space path indicated by the double- 
head arrow in Fig. [Hb). As pointed out in Refi^ the e- 
ph coupling strength in graphene A„k is insensitive to the 
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FIG. 1: (a) Ball-and-stick model of bilayer graphene (Bernal 
stacking). (b) Brillouin zone of graphene and bilayer 
graphene. (c) The electron-phonon coupling strength A„k in 
bilayer graphene versus changing Fermi level -Bp calculated 
along the path (double- head arrow) shown in (b). Solid and 
dashed red lines correspond to A„k of the individual blue and 
red parabolic band in the inset, respectively. The Fermi level 
of neutral bilayer graphene is set at zero, (d) As in (c), for 
each of the two electronic bands of graphite touching at the 
K point (solid blue line) . In (c) and (d) , we show for compar- 
ison the e-ph coupling strength in graphene^^ (indicated by 
the dash-dotted line). 



location of the wavevector k on the Fermi surface. This 
is also the case for bilayer graphene and for graphite. 
Therefore, we drop the index n and the wavevector k 
from now on. In bilayer graphene, the two electronic 
bands near the Dirac point energy exhibit almost identi- 
cal e-ph coupling strengths [Fig.[TJc)]. 

The key factors determining the e-ph coupling strength 
are the density of states around the Dirac point energy 
and the e-ph matrix elements between the initial and the 
final electronic states close to the Fermi level^. The den- 
sity of states of pristine graphene vanishes at the Fermi 
level, whereas bilayer graphene has a finite density of 
states. Despite this difference, at low doping with only 
one band occupied (jii'F — -BdI < 0.2 eV, E-q being the 
energy at the Dirac point and is set to E-q = in the fol- 



lowing discussion), the e-ph coupling strengths in bilayer 
graphene and in graphite are similar (within 5%) to those 
of graphene. This indicates that, as in graphene, there 
is no significant scattering between low-energy electronic 
states in bilayer graphene and in graphite arising from the 
e-ph interaction. This behavior originates from the chiral 
nature of the charge carriers in bilayer graphene^^i^i and 



in graphite 



,25 



it is a matrix-element effect. The dif- 



ference in the e-ph coupling strength between graphene 
and bilayer graphene (or graphite) increases with dop- 
ing. At the largest doping level considered [E-p = 1.5 eV, 
Figs. [T]^c) and HUd)], the couphng strength in bilayer 
graphene and graphite (A — 0.28) is 30% larger than 
in graphene (A = 0.21). As we show in the following, 
these differences result from interlayer interaction. 

In order to determine which phonon modes lead to 
the differences in the e-ph coupling strengths between 
monolayer and bilayer graphene, we decomposed the cou- 
pling strength A„k of both systems into contributions 
from each phonon branch and wavevector. Figures [2jb) 
anddje) show the Fermi surfaces of hole-doped graphene 
and bilayer graphene and the initial wavevector k of the 
electronic state considered. Figures [2{c) and [DJf) show 
the phonon dispersions of both pristine and hole-doped 
graphene and bilayer graphene, respectively. The size 
of the disks superimposed to the phonon dispersions is 
proportional to the contribution to the coupling strength 
A„k arising from the corresponding phonon mode. Fig- 
ures m^c) and Hff) show that both in graphene and in 
bilayer graphene, the major contributions result from 
the highest-energy in-plane vibrations with wavevectors 
connecting the initial and final electronic states on the 
Fermi surface. However, in the case of bilayer graphene, 
the three low-energy optical branches with energy ~ 
10 meV enhance the e-ph coupling strength with respect 
to graphene. The latter vibrations correspond to the 
compression mode (singly degenerate) and to the sliding 
mode of the two layers (doubly-degenerate). 

Since, as in graphene, the e-ph coupling strength in 
bilayer graphene is rather small, even in the heavily 
doped case considered here, it appears unlikely for bi- 
layer graphene to exhibit superconductivity with a tran- 
sition temperature significantly higher than that one may 
expect for graphene. 

Recent angle-resolved photoemission experiments on 
kish graphite^- and on a single crystal of graphitc^^ re- 
ported very different values of the e-ph renormalization, 
namely, A = 0.70 along the KK direction [path indicated 
by the double-head arrow in Fig. [Hb)}^ and A = 0.14 
along the KM direction [Fig. [Hb)}2i. Our calculated e- 
ph coupling strength in undoped graphite (A = 0.034) 
is closer to the estimate of Ref;^. In that work, the 
broadening of the energy distribution curve in the pho- 
toemission spectra was entirely assigned to the e-ph in- 
teraction. This assumption leads to an apparent e-ph 
coupling strength which is enhanced by the contribu- 
tions arising from other interactions. In particular, it 
has been shown that the electronic linewidth in graphite 
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FIG. 2: (a) The electronic energy dispersion and the Fermi level of hole-doped graphene. (b) The Fermi surface (contours) 
and the Brillouin zone (dashed hexagon) of hole-doped graphene. The black dots represent the wavevector k of the electronic 
state considered on the Fermi surface, (c) The phonon dispersion curves of undoped (dashed lines) and hole-doped (solid lines) 
graphene versus the wavevector q along the solid green line shown in (b). The vertical lines indicate the phonon wavevectors 
q such that the final electronic state with wavevector k -I- q is also on the Fermi surface. The size of the disks on top of the 
phonon dispersions is proportional to the contribution of that phonon mode to A„k. (d) to (f): Same quantities as in (a) 
to (c) for hole-doped bilayer graphene but including also interband coupling. The inset of (f) shows one of the three modes 
responsible for the enhancement of the e-ph coupling strength in bilayer graphene. The color (red and blue) and the type 
(solid and dashed) of the curves in (c) and (f) corresponds to the phonon branches in Figs. |4l[a) and|4l[b) andfljc) and|4l^d), 
respectively. 



arising from the electron-electron interaction is sizable^^. 
A similar discrepancy has been pointed out for the case of 
graphene^-'--. However, recent calculations indicate that 
the effect of the electron-electron interaction in graphene 
is not negligible and must be taken into account in the 
analysis of the experimental data^. 

So far we have discussed the effect of the e-ph inter- 
action on the Fermi velocity of the carriers. In what 
follows we focus on the effect of e-ph interaction on the 
phonon linewidths. Figure [Slja) shows the linewidth of 
the doubly-degenerate E2g phonons at the F point and 
of the doubly-degenerate A'l mode at the K point for 
graphene. These phonons exhibit a finite and constant 
linewidth for \E-p\ < Wph/2, where Wph ~ 0.2 eV is the 
optical phonon energy, and a negligible linewidth other- 
wise. The dependence of the phonon linewidth on the 
doping level can be explained by considering that inter- 
band transitions through phonon absorption are forbid- 
den whenever |£^f| > Wph/S. Our calculated linewidth 
of the i?2g phonon is in good agreement with previous 
studies^^i^. 

In bilayer graphene, as a consequence of the inter- 
layer coupling, the four highest-energy modes (originat- 
ing from the F,2g modes of graphene) split into two sets 
of doubly-degenerate Eg and modes^i, with the E^ 



modes 1.1 meV higher in energy than the Eg modes (cf. 
Fig. [3]). Interestingly, at F, only the Eg modes exhibit a 
finite linewidth (1.1 meV) whereas the E„ modes are not 
broadened by the e-ph interaction. It can be shown that 
this difference results from (i) the chiral nature of the low- 
energy electronic state a^^'^^ and (ii) from the fact that 
atoms in the same sublattice but different layers move 
in phase in the Eg modes, while they move out of phase 
in the E„ modest. The calculated linewidths of the Eg 
and the E„ modes are in good agreement with a previ- 
ous study based on a pseudospin effective Hamiltonian 
for the massive Dirac fermions of bilayer graphene^. 

Among the two sets of high-energy zone-center modes, 
only the Eg phonons are Raman active^i. The calcu- 
lated phonon linewidth of the Eg modes can therefore 
be compared directly with the measured broadening of 
the Raman lines. As shown in Fig. [S^b), our calcu- 
lated linewidths are in excellent agreement with those 
reported in a recent experimental study^i. For the pur- 
pose of comparison, we have downshifted the experimen- 
tal linewidths of the E2g phonon in graphene^ and of the 
Eg phonon mode in bilayer graphene^i by 0.6 meV. The 
agreement between our calculations and experiment af- 
ter the subtraction of this uniform background indicates 
that defect-induced scattering and anharmonic effects are 
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FIG. 3: (a) Phonon linewidth in doped graphene for the £2^ 
mode at the F point (solid line) and the A'l mode at the K 
point (dashed line) versus the Fermi level Ep- The filled blue 
squares are the experimental data from Ref.— downshifted by 
0.6 meV (to account for a uniform background), (b) Phonon 
linewidth in bilayer graphene for the Eg mode (solid line) 
and the mode (dash-dotted line) at the F point, and for 
the E mode at the K point (dashed line). The insets show 
one of each of the two doubly-degenerate zone-center modes 
considered here. The filled blue squares are the experimental 
data from Ref.— downshifted by 0.6 meV (to account for a 
uniform background). 



small and similar in magnitude in graphene and bilayer 
graphene. 

We note that the linewidth of the phonons in 
graphene (1.1 meV) and that of the phonons in bilayer 
graphene (1.4 meV) are very similar. This behavior orig- 
inates from the cancellation of the effects of larger elec- 
tron density of states and smaller e-ph matrix elements 
of bilayer graphene as compared to graphene. Because of 
their similar linewidths, the broadening of these modes 
is unlikely to be useful for determining the number of 
graphene layers using Raman spectroscopy. In contrast, 
the linewidth of the highest-energy mode at the K point 
in graphene (the A[ mode) is reduced from 2.9 meV in 
graphene to 1.2 meV in bilayer graphene (the E mode). 
Therefore it should be possible, at least in principle, to 
exploit this difference in two-phonon Raman experiments 
to distinguish between graphene and bilayer graphene. 



Figure Hta) shows the linewidths of the two highest- 
energy phonon branches in pristine graphene. The 
phonon linewidths exhibit maxima at the K point and 
at or near the T point. At the F point, the highest- 
energy phonons decay through electronic transitions with 
no momentum transfer. Off the F point, because of the 
topology of the Dirac cone, non-vertical transitions can 
occur if the wavevector of the phonon is smaller than 
ko = o^ph/fp = 0.035 A {ujph being the phonon energy 
and v-p the Fermi velocity) . These transitions are allowed 
since the phonon wavevectors connect electronic states 
of the same chirality^S. The scattering of phonons with 
wavevector ko = Wph/wp is enhanced because the phase 
velocity of the phonon matches the slope of the Dirac 
cone. Correspondingly, at this wavevector the transverse- 
optical phonon branch exhibits the largest linewidth. 
Unlike the case for the transverse-optical phonons, the 
longitudinal-optical phonons with wavevector ko cannot 
promote electronic transitions, as a consequence of the 
chiral symmetry'^^, and the corresponding linewidth van- 
ishes. 

As shown in Fig. IDJb), in hole-doped graphene, the 
phonon linewidths in the highest energy branches with 
wavevector at the F point or at the K point are negligi- 
ble. However, whenever the phonon wavevector exceeds 
ko in magnitude, intraband electronic transitions can oc- 
cur through phonon absorption. In the case of phonons 
with wavevector close to the K point, this kind of elec- 
tronic transition is suppressed due to chirality^. In hole- 
doped graphene, the Fermi surface consists of two con- 
tours centered around the two inequivalent Dirac points. 
Two maxima are found in the phonon linewidths near K, 
corresponding to the smallest and the largest wavevectors 
connecting electronic states on different contours. 

Figures SJc) and IH^d) show the linewidths of the E^ 
phonons for undoped as well as for hole-doped {Ep — 
— 1.0 eV) bilayer graphene, respectively. In the undoped 
case, the linewidths of transverse phonons with wavevec- 
tor fco are not as large as in graphene since the low-energy 
electronic energy dispersions are nonlinear. In the case 
of hole-doped bilayer graphene, the profile of the phonon 
linewidths is almost identical to the one calculated for 
graphene. This can be explained by considering that 
the electronic density of states per carbon atom near the 
Fermi level in graphene and bilayer graphene are very 
similar. 

The results of our calculations could be confirmed by 
performing detailed inelastic neutron scattering, electron 
energy loss spectroscopy, or inelastic x-ray scattering ex- 
periments. Measurements of this kind have been per- 
formed, for example, in the case of magnesium diboride, 
allowing for a direct comparison with theoretical calcu- 
lations.'*''. 

In summary the Fermi velocity renormalization and 
the phonon line broadening arising from the e-ph in- 
teraction in bilayer graphene and in graphite are stud- 
ied and compared with the corresponding quantities in 
graphene. In bilayer graphene and in graphite, the e- 
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FIG. 4: (a), (b): The phonon linewidth of the highest energy branches in undoped [(a)] and hole-doped [(b)] graphene. The 
color code (red and blue) and the type (solid and dashed) of the line correspond to the phonon branches shown in Fig. EIc) . 
(c), (d): Phonon linewidth of the second highest doubly-degenerate phonon branches at the F point in undoped [(c)] and 
hole-doped [(d)] bilayer graphene. The color code (red and blue) and the type (solid and dashed) of the line corresponds to 
the phonon branches shown in Fig. [2jf). The inset in each panel shows a magnified view of the region near the F point, where 
the symbols represent calculated data points and the lines are a guide to the eye. The vertical dashed line in the inset specifies 
the characteristic wavevector ko (see text). 



ph coupling strength is enhanced by up to 30% at high 
doping as compared to graphene. The calculated doping 
dependence of the phonon linewidth of the zone-center 
Eg mode in bilayer graphene is in excellent agreement 
with recent Raman measurements^. We discussed the 
similarities and the differences in the linewidths of the 
optical phonons in graphene and in bilayer graphene. 
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